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ABSTRACT 

Aims. X-ray Bright Optically Normal Galaxies (XBONGs) constitute a small but not negligible 
fraction of hard X-ray selected sources in recent Chandra and XMM-Newton surveys. Even 
though several possibilities were proposed to explain why a relatively luminous hard X-ray source 
does not leave any significant signature of its presence in terms of optical emission lines, the na- 
ture of XBONGs is still subject of debate. We aim to a better understanding of their nature by 
means of a multiwavelength and morphological analysis of a small sample of these sources. 
Methods. Good-quality photometric near-infrared data (ISAAC/VLT) of four low-redshift (z = 
0.1 - 0.3) XBONGs, selected from the HELLAS2XMM survey, have been used to search for 
the presence of the putative nucleus, applying the surface-brightness decomposition technique 
through the least-squares fitting program GALFIT. 

Results. The surface brightness decomposition allows us to reveal a nuclear point-like source, 
likely to be responsible of the X-ray emission, in two out of the four sources. The results indicate 
that moderate amounts of gas and dust, covering a large solid angle (possibly 4n) at the nuclear 
source, combined with the low nuclear activity, may explain the lack of optical emission lines. 
The third XBONG is associated with an X-ray extended source and no nuclear excess is detected 
in the near infrared at the limits of our observations. The last source is associated to a close (d< 
1 arcsec) double system and the fitting procedure cannot achieve a firm conclusion. 

Key words, galaxies: active - X-rays:galaxies 



* Based on observations made at the European Southern Observatory, Paranal, Chile (ESO Programme ID 
69.A-0554). 
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1. Introduction 



Thanks to the Chandra and 'XMM-Newton surveys, the hard X-ray sky is now probed down to 
a flux limit where the bulk of the X-ray background is almost completely resolved into discrete 
sources (Hasinger et al. 2001 ; Alexander et al. 2003; Bauer et al. 2004; Worsley et al. 2004, 2005). 
Extensive programs of multiwavelength follow-up observations showed that the large majority of 
hard X-ray selected sources are identified with Active Galactic Nuclei (AGN) spanning a broad 
range of redshifts and luminosities. At variance with optically selected quasars. X-ray selected 
AGN are characterized by a much larger spread in their optical properties, especially for what con- 
cerns the intensity of the emission lines. Indeed, a sizable fraction of relatively luminous X-ray 
sources hosting an active nucleus would not have been easily recognized as such on the basis of 
optical observations either because associated with very faint {R > 24) counterparts (e.g., Fiore 
et al. 2003; Mignoli et al. 2004; Civano et al. 2005) or due to the lack of AGN emission lines in 
their optical spectra. The latter class of sources is variously termed as "optically-dull", "optically 
normal" or XBONGs (X-ray Bright Optically Normal Galaxies; Comastri et al. 2002). The com- 
mon meaning of these definitions is that they lack evidence of accretion-driven activity in their 
optical spectra, in contrast with "normal" Seyfert galaxies and quasars. Their X-ray luminosities 

(^ 10"*^ - 10**^ erg s"'). X-ray spectral shape and X-ray-to-optical flux ratio (X/cQ 1) suggest 

AGN activity of moderate strength. Originally discovered in early Einstein observations (Elvis et 
al. 1981) and named optically dull galaxies, the interest on the nature of these sources has gained a 
renewed attention after the discovery of several examples in XMM-Newton and Chandra surveys 
(Fiore et al. 2000; Comastri et al. 2002a,b; Georgantopoulos et al. 2005; Kim et al. 2006). Several 
possibilities were proposed in the literature in order to explain why a relatively luminous, hard 
X-ray source does not leave any significant signature of its presence in the form of emission lines. 

A simple explanation favoured by Moran et al. (2002) and more recently by Caccianiga et al. 
(2007) for faint sources in the Chandra deep fields and brighter object in the XMM-Newton XBS 
survey, respectively, is dilution by the host galaxy starlight. The combination of optical faintness 
and lack of strong emission lines in the observed wavelength range for distant Chandra sources 
or the inadequate observing set-up among brighter nearby XMM-Newton objects (Severgnini et 
al. 2003) may account for the XBONG properties. More in general, if the contrast between the host 
galaxy starlight and nuclear emission is high, AGN emission lines may easily be undetected. The 
physical reason may be ascribed to obscuration, most likely with a large covering factor, or to the 
fact that the lines are not efliciently produced by the central engine. 

If XBONGs are merely obscured AGN, two hypotheses may be envisaged: 

• In order to explain the multiwavelength properties of the XBONG prototype PKS 0312018, also 
known as P3, Comastri et al. (2002) suggested heavy obscuration by Compton-thick gas covering 
almost An at the nuclear X-ray sources. In this way, no ionizing photons can escape to produce the 
narrow emission lines which are observed in "normal" Type 2 narrow-line AGN which are thought 
to have a lower covering fraction following the AGN Unified Scheme (but see Section[3]for a recent 
re-analysis of the X-ray data of P3). 

• According to a detailed multiwavelength analysis of "optically-dull" galaxies in the Chandra 
deep fields, Rigby et al. (2006) conclude that extranuclear dust in the host galaxy plays an 

1 Where X/0 is defined as XjO = log ^ = log/x + C + 
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Alternatively, XBONGs may be members of a class, or classes, of exotic objects for which emis- 
sion lines are either intrisically weak or absent: 

• Radiatively Inefficient Accretion Flows (RIAFs) are expected at accretion rates well below 
those inferred for Seyferts and quasars. A distinctive property of low accretion-rate flows is that 
the standard Shakura-Sunyaev accretion disk is truncated at a relatively large inner radius. As 
a consequence, it cannot generate the "big blue bump" and enough UV photons to photoionize 
the line-emitting circumnuclear gas. The infalling gas is heated to high temperatures and emits a 
hard X-ray power-law by upscattering of low-energy seed photons. According to Yuan & Narayan 
(2004), the SED of source P3 could be reproduced by a RIAF model. 

• XBONGs could be extreme BL Lac objects in which the featureless non-thermal continuum is 
much weaker than the host galaxy starlight. Following Fossati at al. (1998), XBONGs could belong 
to the low-luminosity tail of the blazar spectral sequence based on the anti-correlation between 
luminosity and frequency of the synchrotron peak. 

• A highly speculative possibility is that of a transient AGN phenomenon in the process of 
tidally disrupting a star. If this were the case, the X-ray emission should be witnessing the transient 
accretion phenomenon (see Komossa et al. 2004 for extreme variability events in ROSAT obser- 
vations; see also Gezari et al. 2006 for a luminous flare observed in the GALEX Deep Imaging 
Survey). The transient is most likely over in subsequent follow-up optical observations. 

Finally, it should be noted that diffuse emission from a galaxy group, whose X-ray extended 
emission may have escaped detection in low signal-to-noise X-ray observations, is also possible 
and indeed observed in a few cases (Georgantopoulos et al. 2005). 

While a unique solution may not necessarily hold for all the XBONGs observed in different sur- 
veys, they represent a useful benchmark for a better understanding of the AGN activity and, as such, 
deserve further studies. Ideally, one would need sensitive, high-spatial resolution, multiwavelength 
observations from radio to X-rays. As a first step, in the following we use good-quality photomet- 
ric near-infrared data obtained with ISAAC at VLT of four low-redshift (z = 0.1 - 0.3) XBONGs, 
selected in the HELLAS2XMM survey to search for the presence of a putative nucleus which has 
escaped detection in the optical spectroscopy. The rather obvious advantage of near-infrared data is 
that the effects of dust reddening are minimized and that the nuclear emission in this band should 
rise more rapidly than the stellar light due to the reprocessing by hot dust. At the same time, the 
excellent quality of the near-lR images makes possible to apply a surface brightness decomposi- 
tion technique, already successfully applied by several authors (e.g. Sanchez et al. 2004; Peng et al. 
2006), to search for weak unresolved nuclear emission down to faint near-infrared magnitudes. We 
also discuss the broad-band properties of the four XBONGs using available multiwavelength data. 
Throughout the paper we assume a cosmological model with Hq = 70 km s"^ Mpc"\ = 0.3 and 
Qa = 0.7. 

2. The sample 

In the following, the term XBONG is referred to relatively nearby (z < 0.3), low luminosity (Lx ~ 
2Q42-43 gj.g g-i-)^ X-ray (2-10 keV) selected sources with an optical spectrum typical of early 
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type galaxies. All together they represent a small fraction (~ 5%) of the sources in hard X-ray 
surveys. The four XBONGs were selected from the HELLAS2XMM survey: 3 sources belong to the 
one square degree field (IdF) sample (see Baldi et al. 2002 for a description of the X-ray data and 
Fiore et al. 2003 for the optical identification process), while one object (Abell 1835140) is from 
the survey extension (Cocchia et al. 2007). They represent 40% of the HELLAS2XMM XBONG 
sample (10 objects), and have been selected for infrared observations on the basis of a high quality 
optical spectrum and bright optical magnitudes (R<18). 

The spectroscopic observations were obtained at the ESO 3.6m telescope equipped with 
EFOSC2 during diff'erent observing runs in 2001. The optical spectra, all of them characterized 
by absorption features over a red stellar continuum, are shown in Fig. [1] While the relatively large 
slit width (2 arcsec) used in the spectroscopic observations do not preserve us from a galaxy light 
contamination (and consequent line dilution), the good signal-to-noise ratio (S/Ngy^y^ > 1 3-20) and 
the well suited spectral range (including the expected wavelength of the [OIH]^qq-jx ^^'^ ^'^6563A 
lines, with the exception of Abell 1835140) allow us to place firm constraints on the line fluxes. 

The corresponding luminosity upper limits for the 0[III] line of the four XBONGs (starred sym- 
bols in Fig. |2]l are in the range lO^'-lO"*" erg s about one order of magnitude fainter than that 
measured in the optical spectra of HELLAS2XMM broad-line (BLAGN, filled circles) and narrow- 
line (NLAGN, open circles) AGN with comparable X-ray luminosity (see Fig.|2]and Cocchia et al. 
2007). The upper limits are more than one order of magnitude below the relation between X-ray 
luminosity and 0[III] to X-ray luminosity ratio obtained by Netzer et al. (2006) for a large sample 
of X-ray selected Type 1 and Type 2 AGN spanning the same range of X-ray luminosities. We 
also note that only 0.4% of DR6 SDSS AGN (about 7000 quasars; Adelman-McCarthy et al. 2007) 
have 0[III] luminosity below 10""' erg s 

The redshifts obtained from the position of the principal absorption lines are 
z = 0.319, 0.159, 0.154, 0.251 for PKS 03120017, PKS 03120018, Abefl 2690013 and 
Abell 1835140, respectively (Fiore et al. 2003; Cocchia et al. 2007). The X-ray unabsorbed 
luminosities are reported in Table [1] Radio observations at 5 GHz are available for the two objects 
in the PKS 0312-77 field. Source PKS 03120017 was detected with a flux density of 1.3 mJy at 

5 GHz, while source PKS 03120018 remains undetected (Fschz < 015 mJy at the 3cr level). An 
observations at 8.7 GHz is also available for source PKS 03120017 which was detected with a flux 
density of 0.5 mJy. 

3. X-ray spectral analysis 

The results of the spectral analysis of the X-ray sources in the IdF sample are presented in Perola 
et al. (2004) with particular emphasis towards the statistical properties of the sample. The X-ray 
spectra of the four XBONGs discussed here were re-analyzed including an additional observation of 
the Abell 1835140 field (not included in the Perola et al. paper) using the most recent calibrations. 
The data were processed using standard sa^^ v6.1.0 and ftools tasks. 

The pn spectra were grouped with a minimum of 10 counts per bin (for PKS 03120018 and 
PKS 03120017) and 15 counts per bin (for AbeU 26900013 and AbeU 1835140) using the task 
GRPPHA, and fitted with xspec v 11.3.0 (Arnaud 1996), while MOSl and MOS2 data, having much 



^ XMM-A'evv?o;i Science Analysis System; see |http://xnim.vilspa.esa.es/extemal/xnim_sw-cal/sas-frame.shtml| 
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Fig. 1. Rest-frame optical spectra of the sources PKS 03120017 and PKS 03120018, Abell 2690013 
and Abell 1835140, taken with the ESO 3.6m telescope. The stellar absorption lines are labeled. 
The Hot line is also labeled at the expected wavelength. 
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Fig. 2. The [OIir\^QQ-j^ luminosity versus the 2-10 keV luminosity for 59 sources of the 
HELLAS2XMM sample as presented in Cocchia et al. (2007). Starred symbols mark the four 
XBONGs discussed in this paper; filled circles - broad-line AGN; open circles - narrow-line AGN; 
open triangles and open squares - emission-line and early-type galaxies, respectively. 

lower counting statistics, were left unbinned and fitted with the C-statistic (Cash 1979). We 
checked and found consistency between the pn and MOS data. However, the X-ray spectral results 
reported in Table [T] and the best-fitting spectra shown in Fig. [3] are referred to the pn data alone, 
given the lower counting statistics of the MOS data. In the following we summarize the main 
X-ray spectral results. 

Two sources, PKS 03120018 and Abell 1835140, are best fitted by an absorbed power law with 
photon index in the range F a 1 .5-2. 1 and column densities (Nh ~ 4-7x10^' cm"^), typical of a 
mildly obscured AGN. Despite the relatively large errors in the determination of the X-ray spectral 
parameters, this analysis seems to rule out Compton-thick absorption, at least in the energy range 
covered by XMM-Newton. The radio-loud source PKS 03120017 is best fitted by a relatively steep 
F ^ 2.4 spectrum (but consistent with F ~ 2) and no absorption in excess to the Galactic value (see 

Table 1. X-ray spectral analysis parameters of the pn data. The X-ray luminosities are corrected 
for absorption. For source AbeU 1 835 140, the spectral parameters of the two available observations 
are reported. 



Src. ID 


kT 


r 


Nh 




F2-IO ke 


V 


^2-10 keV 




(keV) 




(cm"-) 




(erg cvcT^ 


s-') 


(erg s-') 


PKS 03120017 




7 4.+0.6 


< 2.1 X 10-' 


3.7/11 


1.4x10" 


14 


5.0xlO'*2 


PKS 03120018 




9 1+0.8 
^•1_0.4 


e.rll X 10^1 


9.6/9 


1.5x10- 


14 


l.lxlO'*^ 


Abell 26900013 


0.73±0.06 






34.4/25 


1.5x10- 


14 


9.6x10'*' 


Abell 1835140 






4.3!|5 X 10" 


14.0/12 


5.3x10- 


14 


9.0x10''^ 






^•'-0.3 


6.5!|^ X 1021 


32.9/35 


1.8x10- 


13 


3.3xlO« 
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Fig.3. XMM-Newton (pn) spectra of the sources PKS 03120017, PKS 03120018, Abell 2690013 
and Abell 1835140 (from top-left to bottom-right), fitted in the 0.5-8 keV range. The solid line 
represents the best-fitting model. The spectral parameters are reported in Table [T] For source 
Abell 1835140, the spectra of the two available observations are reported. 

Table[B. 

Abell 26900013 is well fitted using a thermal (mekal) model at soft energies and a weakly 
constrained power-law tail at higher (>2 keV) energies. The thermal spectrum is consistent with 
the extended appearance of the X-ray emission, as shown by the overlay of the X-ray contours on 
the Kj-band image (see Fig.]?]). 

Significant X-ray flux variability (a factor about 3) is present between the two XMM-Newton 
observations of Abell 1835140. This object was also observed twice and detected by Chandra (in 
1999 and 2000, for « 20 and 10 ks, respectively). The Chandra data (especially the data of the first 
and longer observation) confirmed the XMM-Newton spectral results and allowed for a long-term 
variability study of this source. Over ^ 3 and half years, Abell 1835140 was found to be variable 
by a factor of w 7.5, with its 2-10 keV flux increasing from the value of ^ 2.4 x 10"'"* erg cm"^ s"^ 
in 1999 to * 1.8 x 10 erg cm"^ s ' in 2003, with no evidence for significant changes in the X-ray 
spectral parameters. Even if the flux variability could be connected with the tidal disruption of a 
star, the increasing of the X-ray flux allows us to rule out such a possibility. 

4. Near-Infrared data analysis 

4.1. Near-Infrared imaging observations 

Deep near-infrared observations have been carried out using the Infrared Spectrometer and Array 
Camera (ISAAC) mounted on the ESO VLT (Very Large Telescope). ISAAC is equipped with a 
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Fig. 4. X-ray contours obtained from the 0.5-10 keV adaptively smoothed pn image of source 
Abell 2690013 overlaid on the K^-band image. The image is 60"x60". 



Rockwell Hawaii 1024x1024 HgCdTe array, with a pixel scale of O'.'147/pixel and a field of view 
of 2^5 X 2^5. The observations have been collected in service mode over several nights during 
September 2002 undergood seeing conditions (<0'.'8). The ISAAC images were taken in two near- 
IR bands (the J and Kjjfilters) and were reduced using the DIMSU1V0 package, following standard 
procedures. A detailed discussion on the data reduction of similar images taken for other sources 
revealed in the HELLAS2XMM IdF is reported in Mignoli et al. (2004). The total exposure times 
for each galaxy is 600 sec in J and 1800 sec in K, band. The effective seeing on final frames 
ranges between 0.5 and 0.8 arcsec (FWHM). The magnitude limits for point-like sources are about 
/(Um)=22 and Kj(lim)=21.5 in all the images. The Kj-band images of all the four sources are 
shown in Figure |5] 



4.2. The fitting procedure 

In order to search for the presence of a putative nucleus, we applied the surface brightness decom- 
position technique using GALFIT (Version 2.0.3b, Peng et al. 2002), a two-dimensional algorithm 
designed to extract structural parameters directly from galaxy images combining several analytical 
models and by convolving them with the image Point Spread Function (PSF). The PSF has been 
modeled adopting several field stars in each ISAAC image. The fitting has been performed in a box 
around each object large enough (220x220 pixels) to contain most of the bulge. 

The fitting procedure described below was performed separately in both J and K, bands, al- 
lowing us to cross-check the results and test for a possible color trend. In order to better constrain 
the galaxy structural parameters (magnitude and effective radius), first we fitted each object in the 
J band, where the AGN contribution to the total light was expected to be lower (according to a 

^ The Kj (K-short) filter is centered at shorter wavelenght than the standard K filter in order to reduce 
thermal background. 

Deep Infrared Mosaicing Software, developed by P.Eisenhardt, M.Dickinson, A.Stanford and J.Ward., 
and available at the site ft p : // iraf . noao. edu/ contr ib/ dimsumV2 , 
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Fig. 5. K,-band images: PKS 03120017, PKS 03120018, Abell 2690013 and Abell 1835140, from 
top left to bottom right; the size of the images is 20"x20", 15"xl5", 30"x30", 15"xl5", respec- 
tively. 



standard AGN SED), as actually found in our analysis. Then we appUed GALFIT to the Kj-band 
images, using the galaxy parameters obtained in the J band as start-up values for the nuclear de- 
composition. 

As a first guess of the profile of the host galaxy, we applied a Sersic model (Sersic 1968) of the 
form: 

Mr)=)U.e-4(^)""-i] (1) 

where /ig is the effective surface brightness, the effective radius, n the Sersic index, and k a 
constant determined from n in order to have half of the integrated flux within r,,. The Sersic index 
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assumes the value n ~ lin case of exponential profile and value n = 4 in case of a de Vaucouleurs 
profile (de Vaucouleurs 1948). A fit with a Sersic model returns n > 6 in all the sources. Since this 
index is higher than the typical values found in early-type galaxies (e.g. Andredakis et al. 1995; 
Cassata et al. 2005) and given the early type like galaxy spectra of our sources, we adopted a 
de Vaucouleurs profile, hereafter referred to as MODEL 1 . 

If residual emission in the innermost region were still present after fitting the galaxy profile 
with MODEL 1, we added the contribution of a point-like source (MODEL 2). The latter is mod- 
eled using the PSF profile obtained by averaging stars in the field. The centroid of the unresolved 
component and that of the host galaxy have been linked each other in the fitting procedure. We 
considered the unresolved component as detected if (1) the;^^^ significantly improves with the in- 
clusion of the nuclear component (see Section |43T l and (2) it is brighter than the image magnitude 
limit for point-like sources. The results of the fitting are shown in Table|2] 

In order to check the GALFIT capabilities to detect a faint unresolved source at the center of a 
bright galaxy, accurate simulations were performed using the IRAF task mkobject. We started with 
a r'^"^ profile and then a point-like source has been added. In the simulations, we adopted a grid 
of magnitudes in the range Kj=17 - 22 for the point-like component and Kj=14 - 22 for the host 
galaxy; different effective radii (from 30 to 120 pixels) for the host component have been used. 
The input parameters have been chosen in an appropriate range to reproduce the observed data and 
to span a wider but realistic parameter space. The simulated images have been convolved with the 
seeing and the noise of our images; the PSF has been modeled with the parameters determined from 
the field stars. We applied GALFIT on the simulated sources first modeling with a de Vaucouleurs 
profile and then adding a central unresolved component, as performed for the XBONGs. The fitting 
results are fully consistent with the input parameters. Uncertainties (A mag ~ 0.02) on the determi- 
nation of the nuclear magnitude have been found in the case of a very faint (K, >22.5) unresolved 
source hosted by a bright (K, ^15) galaxy. We conclude that at the magnitude hmit of our images 
for point-like sources, the errors on the parameters due to the GALFIT detection capabilities are 
not affecting our results significantly. 



4.3. The fitting resuits 

The best-fitting parameters (m/,oj,, radius, m„uck'us) obtained in each band for MODEL 1 and 
MODEL 2 are reported in Table |2] The fifth column of Table |2] gives the effective radius, which 
is always larger than the seeing disk (up to 3.5 kpc at z = 0.3). Since the uncertainties quoted by 
GALFIT are very small (see Peng et al. 2002), the errors reported in Table|2]have been estimated 
comparing the best-fitting results obtained using the different stars in the images one by one as a 
model for the PSF; this method has been applied because it is known that steep profiles are very 
sensitive to the accuracy of the used PSF (see the online manua|fl). 

• In two out of the four sources (PKS 03120017 and PKS 03120018), we found clear evidence 
for an unresolved component on top of the elliptical galaxy profile. The Kj-band residual images, 
obtained by subtracting the fitted models from the original images, are shown in Fig. |6j the resid- 
uals obtained with MODEL 1 are shown in the left panels, while those obtained with MODEL 2 
are shown in the right panels. Introducing an unresolved component improves the quality of the fit, 

^ http://zwicky.as.arizona.edu/ cyp/work/galfit/TFAQ.html. 
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Table 2. Results of the fits performed with GALFIT in the J and in the K, band; the two entries for 
each source in each band correspond to the fitting with MODEL 1 and with MODEL 2, respectively. 



id 


band 


Model 


'-'-'-no SI 








d.o.f. 


y- /d.o.f 


PKS03 120017 


J 


1 


16 94+0 03 


4.79+0.51 




45983 


48835 


0.94 




J 


2 


16 92+0 05 


5 64+0 86 


20 75+0 29 


45614 


48832 


0.93 




K 


1 


15.24+0.01 


4.96+0.17 




54680 


48835 


1.12 




K, 


2 


15.21+0.01 


6.50+0.10 


18.39+0.24 


49722 


48832 


1.02 


PKS03120018 


J 


1 


16.51+0.02 


3.73+0.30 




75401 


48835 


1.54 




J 


2 


16.52+0.01 


4.03+0.40 


20.30+0.35 


67046 


48432 


1.37 




K, 


1 


15.25+0.01 


3.30+0.05 




85728 


48835 


1.75 




K, 


2 


15.20+0.01 


5.04+0.20 


18.09+0.10 


71944 


48832 


1.47 


Abell2690013 


J 


1 


15.41+0.01 


7.34+0.10 




51227 


48835 


1.05 




K, 


1 


14.29+0.01 


6.86+0.20 




106194 


48835 


2.17 



mostly in the Kj band, as evident from the residual images and the^^ values reported in Table |2] 
The excess with respect to the bulge profile in the innermost region increases towards longer wave- 
lengths. The variation in the reduced obtained by adding a nuclear component, is significant at 
a high (> 99.9999%) confidence level, estimated using an F-test, for both sources. We note that 
the statistic is computed over a very large number of degrees of freedom (see Table |2]l, corre- 
sponding to each individual pixel in the image. In order to further assess the improvement in the;^'^ 
statistic obtained by the addition of a nuclear component, whose contribution is relevant only over 
a small portion of the image, we ran GALFIT (both MODEL 1 and MODEL 2) over the central 
innermost part of the image (25x25 pixels) keeping fixed the best fit parameters obtained by the 
analysis of the entire image. Also in this case, the fit quality improves significantly (> 99.9999% 
level, according to the F-test) using MODEL 2 with respect to MODEL 1, indeed enforcing the 
findings obtained considering a larger region. 

• The fitting procedure found no evidence for a unresolved nuclear component in 
Abell 2690013. The residual images obtained with MODEL 1 in both J and Kj bands show no 
need for a central unresolved source in the fitting model; the residuals in the band are shown in 
Fig El If a nuclear component is present, its magnitude is fainter than the limit of our observations. 
The absence of a central unresolved component is also supported by the X-ray analysis: the X-ray 
emission is well fitted by a thermal model consistent with the extended appearance, as shown in 
Figure m 

• Thanks to the high quaUty of the ISAAC images with sub-arcsec seeing, it has been possi- 
ble to unveil the complex morphology of source Abell 1835140. It is a close interacting system, 
composed by two nearby bulges embedded in a common envelope (see bottom right panel of Fig. 
|5]l. The angular separation of the two nuclei is 0.9 arcsec. The morphological fitting was attempted 
including simultaneously the two components, despite the small angular separation. We fitted a 
Sersic model for both galaxies: in the J band the fit did not provide a statistically reasonable so- 
lution (i.e., by means of^^), while in the K, band we obtained a steep profile for the northern 
nucleus and a flatter profile for the southern one. Therefore we performed a new fit, adding an 
unresolved component only to the northern source and a de Vaucouleurs profile to both the bulges. 





Fig. 6. Residual images in the K,- band (galaxy - model) of PKS 03120017 (top) and 
PKS 03120018 (bottom) obtained applying a model without (on the left) and with (on the right) a 
central unresolved component. 
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Fig. 7. Residual image (galaxy - model) in the K,- band of source Abell2690013 obtained applying 
MODEL 1. 
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Fig. 8. Near-infrared color obtained for the AGNs hosted in PKS 03120017 and PKS 03120018 
compared with literature data. Filled (red) dots are the optically selected quasars with near-infrared 
counterparts in the 2MASS (Barkhouse et al. 2001), while the open (black) squares are the 2MASS 
selected red quasars from Hutchings et al. (2003). The tracks plotted represent the color-redshift 
relations for a composite quasar template with different extinction values {E(B - V)-Q, 0.5, 0.8, 
from bottom to top). 

as performed for the other XBONGs. The results are suggestive of the presence of a point-like com- 
ponent in the northern source of the double system. Since the fitting did not provide a statistically 
acceptable description of the two sources, leaving significant residuals because of their low spatial 
separation, this result should be considered just qualitative. 

4.4. Near-Infrared color 

The surface brightness decomposition technique allowed us to separate the galaxy from the unre- 
solved nuclear contribution in two out of four XBONGs. The host galaxy parameters (m/,£,s,, r^,) and 
the 7-Kj colors are in good agreement with the Kj-band luminosity vs. radius relation (Pahre 1999) 
and the color typical of elliptical galaxies (Cutri et al. 2000; Mannucci et al. 2001), respectively, 
lending further support to the morphological analysis. The total luminosity of the galaxies in our 
sample is of the order of L* of early-type galaxies in the band (Kochanek et al. 2001). 

The nuclear near-infrared colors of PKS 03120018 (7-K,=2.21 ± 0.36) and PKS 03120017 
(7-K5=2.36 ± 0.37) are reported in Figure [8] along with Uterature data. Dots represent the optically 
selected quasars with near-infrared counterparts in the 2MASS (Barkhouse et al. 2001), while the 
open squares are the 2MASS selected red quasars from Hutchings et al. (2003) in a redshift range 
(0.1 < z < 0.3) similar to that of our sources. Even if affected by large uncertainties, mainly due 
to the extremely faint 7-band magnitudes, the J-K, colors of the two XBONGs are consistent with 
those of 2MASS red quasars at the same redshift. The tracks plotted in Fig. [8] represent the color- 
redshift relations for a composite quasar template from the Large Bright Quasar Survey (LBQS; 
Francis et al. 1991) extended to the near-lR using the mean radio-quiet quasar energy distribution 
(SED) by Elvis et al. 1994 (see also Maiohno et al. 2000). We have included the eff'ect of dust 
attenuation on the LBQS template using a dust screen model and the SMC extinction law (Pei at al. 
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1992), to keep into account the presence of nuclear or galactic dust. The curves represent different 
extinction values iE(B - V)-0, 0.5, 0.8 from bottom to top). At the face value, the best-fit J-Ks 
color of the nuclei implies E{B -V)- 0.5 - 0.8 which corresponds to Nh - 3 - 5 x 10^' cm"^ if 
we use a standard Galactic dust-to-gas ratic^ as a conversion factor. 



5. Simulated optical spectra 

In order to check whether the extinction values, estimated from the nuclear near-infrared colors (J- 
Ks), could explain the lack of the optical emission lines, we compared the observed optical spectra 
of our XBONGs with simulated optical spectra. We adopted the SDSS quasar composite spectrum 
(Vanden Berk et al. 2001) extended to the near infrared (see previous section) to reproduce the 
nuclear component. 

The approach we used is as follows: 

1 . using the IRAF task de redden, the quasar template was reddened with different extinction val- 
ues: EiB - y)=0, 0.2, 0.4, 0.6, 0.8; 

2. the templates were redshifted to the redshift of the sources using the IRAF task dopcor; 

3. the spectra were normalized to the K, band on the basis of the estimated magnitude for the 
unresolved central source; 

4. the nuclear simulated spectra were summed to the observed XBONG optical spectra. 

We computed simulated spectra for both PKS 03120017 and PKS 03120018. In spite of the good 
signal-to-noise ratio of the observed spectra, the [O/Z/Jjooy^ emission line is not visible in all the 
summed spectra, even in the case without reddening. On the contrary, it is possible to reveal the 
presence of an emission line at the expected wavelength of the redshifted Ha^^^^^^k '^^e spectra 
with E{B - y)=0. 

In order to quantify the emission-line detection limits as a function of E(B - V), we measured 
the equivalent width (EW) of the Ha in the simulated spectra. The average EW of the Ha line in 
the SDSS QSO composite spectrum is ~ 190 A (Vanden Berk et al. 2001). Since the EW is defined 
as the ratio between the line and the continuum fluxes, its value does not change applying the ex- 
tinction in the nuclear simulated spectra. On the contrary, the EW in the total spectra depends from 
both the nuclear and the steady host galaxy continuum, thus it decreases for increasing extinction. 
Figure|9]shows the equivalent width measured in the summed spectra as a function of the E(B - V) 
along with the 3cr level of detectability (dot-dashed line) for both sources. The EW detection limit 
has been estimated following the relation presented in Mignoli et al. (2005): 
SLxA 

EW{detection limit) - — (2) 

(S /N)cont 

where A is the resolution element (in A) of our spectra, SLis the significance level of the detectable 
line, expressed in terms of cr above the continuum noise, and S /N is the signal-to-noise ratio of the 
continuum in the total spectra. The extinction values for which the Ha equivalent width is lower 
than the estimated detection limit is E(B-V) ~ 0.5, in fair agreement with the values found with the 
near-infrared morphological analysis. In the next section we will compare all the extinction values 
found so far. The above described results are not strongly dependent on the assumed Ha EW. 

^ £(B - V) = 1.7 X 10-^^ X Nh cm^^ (Bohlin et al. 1978). 
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Fig. 9. The observed equivalent width of the Ha line plotted against the E{B - V) of the sim- 
ulated spectra; filled squares in the upper panel represent PKS 03120018 and stars in the lower 
panel PKS 03120017. The dot-dashed lines represent the EW detection limit curves as estimated 
following Mignoh et al. (2005). 

In fact, also adopting a lower value (70 A, corresponding to the 2cr lower bound of the average 
intensity of FBQS survey; Brotherton et al. 20010), a weak line is detected in the simulated spectra 
with null extinction. 

6. Discussion 

The surface brightness decomposition technique allowed us to detect and separate the nuclear from 
the galaxy contribution in two (PKS 03120018 and PKS 03120017) out of the four XBONGs. 

There is no evidence for the presence of a significant nuclear component in Abell 2690013; the 
magnitude and the effective radius of the galaxy are in good agreement with the Kormendy relation 
extended to the K^ band (Pahre 1999). Moreover, the presence of an active nucleus could be ruled 
out also thanks to a more careful analysis of the X-ray image. The X-ray source is clearly extended 
(Fig. nil and the X-ray spectrum is well fitted by a thermal model. 

The presence of two close {d <0.9 arcsec) galaxies embedded in a common halo with disturbed 
morphology in source Abell 1 835 140 prevented us from a detailed analysis with GALFIT. A qual- 
itative analysis suggested the presence of a unresolved component in the northern source of the 
system, possibly associated with the X-ray emission. 

The estimated nuclear 7-Ks colors, even if affected by large errors due to the uncertainties 
on the nuclear 7-band magnitudes, are consistent with extinction values of E(B - V) ~ 0.5 and 
0.8 for PKS 03120018 and PKS 03120017, respectively. For both sources, the lower limit found 
with the optical simulated spectra, as explained in Section |5] is E{B - V) ~ 0.5. These values 
are in agreement with the gas absorption column densities found with the X-ray spectral analysis, 
assuming a Galactic dust-to-gas ratio (Bohlin et al. 1978). 

The SEDs of PKS 03 120017 and PKS 03 120018 are reported in Figure[TO]and[n] Open squares 
represent the host galaxy and filled squares the nuclear emission. The X-ray spectra are the same as 
^ The dispersion around the average Ha value is not quoted for the SDSS composite spectrum. 
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Figure |3] taking into account the errors on the X-ray flux. We note that source PKS 03120017 has 
been detected at 5 and 8.7 GHz, with flux densities of 1.29 and 0.55 mJy, respectively. In the left 
panels the data are compared with the Type 1 AGN template from Elvis et al. (1994), normaUzed to 
match the 5 keV rest-frame luminosity for both sources as well as the 5 GHz rest frame luminosity 
for PKS 03120017. The X-ray to near-infrared ratio of both sources is consistent with that of a 
Type 1 AGN. The radio spectral index of source PKS 03120017 is very steep (a ~ 1.7), at variance 
with the value expected for an active nucleus (a < 0.8). However, the radio observations are not 
simultaneous and thus source variability may afifect the slope estimated above. 

The detection of PKS 03 120017 in the radio band suggested (Brusa et al. 2003) that it could be 
an example of a high-energy-peaked BL Lac (the so called HBL). To shed light on the nature of 
the broad-band emission of this object, we report in Fig. [10] (right panel) the SED with the lowest 
radio luminosity belonging to the "blazar sequence" (Fossati et al. 1998), roughly matching the 
observed 5 GHz flux density of our source. The X-ray-to-radio luminosity ratio of PKS 03120017 
is one order of magnitude lower than that expected for a HBL. The HBL SED does not adequately 
fit the nuclear data of PKS 03 120017 suggesting that either the blazar sequence breaks down at low 
luminosities or, more likely, that PKS 03120017 is not a BL Lac object. 

The presence of Compton-thick absorption for source PKS 03120018 was favoured by 
Comastri et al. (2002), on the basis of the comparison of the source SED (including the host galaxy) 
with that of NGC 6240, the prototype of this class of objects (see their Fig. 5). The near-infrared 
nuclear fluxes, obtained with the present analysis, along with a re-analysis of the X-ray spectrum, 
consistent with Compton-thin obscuration (Nh~5 xlO^' cm"^), make the Compton-thick hypothe- 
sis no more tenable (Fig.[TTl right panel). 

A multiwavelength spectral fit to the total emission of source PKS 03120018 with a RIAF plus 
thin-disk model was presented by Yuan & Narayan (2004). Even if a detailed analysis of the RIAF 
model is beyond the scope of this paper, their model is also reported in Fig. [TT] (right panel). We 
note that the solution proposed by Yuan & Narayan (2004) does not provide a good fit to the nuclear 
data; however, given the large number of free parameters in RIAF models, we cannot rule out that 
inefficient accretion may explain the PKS 03120018 SED. 

To explain the properties of more than half of the XBONGs selected in the Chandra deep fields 
at high redshift, Rigby et al. (2006) proposed the presence of extranuclear dust and gas distributed 
on large scales (i.e., kpc dust lanes). Dust lanes have been observed in HST images of nearby 
galaxies hosting a Seyfert 2 nucleus (Malkan et al. 1998) and seem also to be ubiquitous in nearby 
early-type active galaxies, selected from the Ho et al. (1995) sample and analyzed by Lopes et 
al. (2007). X-ray observations of galaxies belonging to the Malkan et al. sample suggest a strong 
correlation between the presence of dust lanes and X-ray absorption (Guainazzi et al. 2001, 2005). 
Moreover, the presence of extranuclear dust lanes in both samples is always associated with optical 
emission lines, at variance with the XBONG definition. If extranuclear dust were responsible of the 
XBONG classification, as proposed by Rigby et al. (2006), it should be visible in our sources in the 
optical band, since our XBONGs are much brighter and at lower redshift. We have retrieved from 
tiie HST archive two ACS/HRC images of source PKS 03120018 taken in the F475W and F625W 
filters. The optical images along with near-infrared images do not show any obvious large-scale 
dust structure. 
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Fig. 10. The broad-band spectral energy distribution of source PKS 03120017 (filled 
squares=nuclear data; open squares=host galaxy data) compared with the Type 1 AGN SED (Elvis 
et al. 1994; left panel) and with a blazar SED with radio luminosity consistent with that of our 
object (Fossati et al. 1998; right panel). The AGN SED is normalized both to the radio and the 
X-ray data (rest frame 5 GHz and 5 keV, respectively). 
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Fig. 11. The broad-band spectral energy distribution of source PKS 03120018 (filled 
squares=nuclear data; open squares=host galaxy data) is compared with the Type 1 AGN SED 
(Elvis et al. 1994) normalized to the X-ray data (5 keV, left panel) and with that of the highly ob- 
scured Seyfert 2 galaxy NGC 6240 (continuous line) and with the RIAF model (dot-dashed line; 
Yuan and Narayan 2004, right panel). 

If heavy nuclear absorption or absorption related to the presence of dust lanes are not respon- 
sible for the lack of optical emission lines in PKS 03120017 and PKS 03120018, our results re- 
quire an alternative explanation. One speculative possibility is that the ionizing flux does not reach 
the narrow line regions because absorbed by intervening dust and gas, with E{B - V) - 0.5 and 
Nh ~10^'cm"^, covering An at the central engine. 

Having estimated the host galaxy magnitude, it is possible to compute the black hole mass 
adopting the Marconi & Hunt (2003) relation, safely assuming that it holds at the redshift of our 
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Fig. 12. Black hole masses (in units of solar mass) versus 2-10 keV luminosity. Filled squares are 
the two XBONGs, open squares the LLAGN from Panessa et al. (2006). The lines represent the ex- 
pected correlation between the two plotted quantities for two different assumptions on the bolomet- 
ric correction (k\-,oi - 15 and kboi - 30, solid and dashed line, respectively) and L/L^^/^^O. 001-0.1 
(from upper left to bottom right). 



sources. The rest frame K^-band luminosity has been computed using an evolving elliptical galaxy 
template (Bruzual & Chariot 2003) to properly account for the K-correction. The black hole 
masses obtained are 6.5 x 10**Mo and 1.4 x IO'^Mq for PKS 03120017 and PKS 03120018, respec- 
tively. Assuming a range of the bolometric correction value (15-30) which, following Marconi et 
al. (2004; see also Elvis et al. 1994), is appropriate for the X-ray luminosity of the XBONG nuclei 
(5.0 X 10^*2 erg sec"' and 1.1 x 10'*^ erg sec"' for PKS 03120017 and PKS 03120018, respectively), 
it is possible to derive the bolometric luminosity and, therefore, to estimate the Eddington ratio 
{LboilLEdd) of the two sources. The resulting Eddington ratio is close to LhoilLEdd ~ 10"^ for both 
sources. Should the bolometric correction be significantly lower (as it seems to be the case among 
low luminosity AGN e.g. Ho 1999) then the estimated Eddington ratio would correspondingly 
decreases. At such low values of the Eddington rate, the accretion is expected to be radiatively 
inefficient and described by a RIAF solution. However, it has also been suggested that the thin disk 
solution is stable down to Eddington ratios as low as 10"^ (Park & Ostriker 2001). Therefore, the 
lack of optical emission lines may be explained by a either a RIAF model or by standard accretion 
disk which, at low values of the Eddington ratio, does not produce enough UV photons able to 
ionize the narrow line regions. 

In Fig.[12]the black hole masses of the two XBONGs as a function of the X-ray luminosity (filled 
squares) are reported along with the low-luminosity AGN (LLAGN) from Panessa et al. (2006). 
Keeping into account the uncertainties on our estimates, the black holes that reside in our objects 

^ The galaxy template includes an exponential star formation with a time-scale t=0.3 and formation red- 
shift ~6. 
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are more massive than those of the LLAGN of Panessa et al. (2006) and, for the same X-ray 
luminosity, are characterized by a lower Eddington ratio. 

7. Conclusions 

We have presented a multiwavelength analysis of four XBONGs selected from the HELLAS2XMM 
survey. For these sources, deep near-infrared images taken with ISAAC at VLT, good-quality op- 
tical spectra and XMM-Newton data are available. Applying the morphological decomposition 
technique, we were able to detect the presence of a nuclear component in two out of the four 
sources (PKS 03 1 200 1 7 and PKS 031200 18). There is no evidence of nuclear emission in the near- 
infrared in Abell 2690013; moreover, the X-ray appearance is consistent with an extended source. 
For source Abell 1835140, the near-infrared images reveal a complex morphology, where two 
sources are embedded in a common envelope. The main issue about the XBONG nature is whether 
they represent a truly distinct class or, rather, they are a mixed source population. Our results point 
towards the latter hypothesis. 

The results regarding the nature of the XBONGs with nuclear component and the lack of optical 
emission lines can be summarized as follows. 

- Source PKS 03 1 200 1 7 and PKS 03 1 200 1 8 are well described by a mildly obscured (E{B - V ) = 
0.5 - 0.8) optically weak nucleus, responsible for the X-ray emission, hosted by a bright galaxy 
(mag™'''-mag'£" ~ 4). 

- The lack of optical emission lines cannot be attributed to observational limitations such as a 
inadequate observational setup or low signal-to-noise-ratio, at least for the objects in the present 
sample for which high-quality spectroscopic and photometric observations are available. 

- We can safely discard for the two objects both a Compton-thick scenario (as found by 
Caccianiga et al. 2007 for a sample of elusive AGN in the XBSS) and an important blazar 
contribution. A RIAF solution seems to be supported by the estimated values of the Eddington 
ratios in the two objects and is not ruled out by the broad-band SED fitting, also because of 
the large number of free parameters in the RIAF model that can be tuned in order to reproduce 
the observed SEDs. A weak nuclear source, described by a standard accretion disk solution, 
but not powerful enough in the production of UV photons, would also provide an acceptable 
description of the observations. 

- The presence of a thin nuclear gas and dust structure (as argued by Cocchia et al. 2007) 
covering 4jt at the nuclear source, combined with the low level of activity of the BH, could 
prevent the ionization of the narrow-line regions and produce also the extinction we measured. 

The analysis of already obtained VLT/VIMOS-IFU (Integral Field Unit) observations could 
help in the detection of the typical AGN emission lines down to a very faint flux limit and in 
the study of the gas kinematics. 
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